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j>! Studies of the diffuse X-ray emitting gas in galaxy clusters have provided pow- 

erful constraints on cosmological parameters and insights into plasma astro- 
physics. However, measurements of the faint cluster outskirts have become 
possible only recently. Using data from the Suzaku X-ray telescope, we deter- 
mined an accurate, spatially resolved census of the gas, metals, and dark mat- 
ter out to the edge of the Perseus Cluster. Contrary to previous results, our 
measurements of the cluster baryon fraction are consistent with the expected 
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universal value at half of the virial radius. The apparent baryon fraction ex- 
ceeds the cosmic mean at larger radii, suggesting a clumpy distribution of the 
gas, which is important for understanding the ongoing growth of clusters from 
the surrounding cosmic web. 

Galaxy clusters provide critical constraints on cosmological parameters which are indepen- 
dent from those determined using type la supernovae, galaxy surveys, and the primordial cosmic 
microwave background radiation (CMB) (1-3). In particular, knowledge of their baryon content 
is a key ingredient in the use of clusters as cosmological probes (4, 5). Most baryons reside in 
the hot, diffuse, X-ray emitting intra-cluster medium (ICM). Until recently, X-ray observations 
have generally targeted only the inner parts of clusters, where the emission is brightest, leaving 
a large fraction of their volumes practically unexplored. Estimates of the gas and total masses 
at large radii have relied on simple model extrapolations of the thermodynamic properties mea- 
sured at smaller radii. 

X-ray spectroscopy of the outer regions of galaxy clusters was made possible only recently 
using the Suzaku satellite. Because of its much lower instrumental background than other X- 
ray observatories, Suzaku can measure the characteristics of the faint emission from cluster 
outskirts more reliably. Even so, few such observations have been published, and the thermo- 
dynamic profiles at large radii are not well resolved (6-11). The Perseus Cluster of galaxies 
is the brightest, extragalactic extended X-ray source. It is a relaxed system, both closer (at a 
redshift of 0.0183) and substantially brighter than any of the other clusters for which Suzaku 
has been previously used to study the ICM properties at large radii. Importantly, its large an- 
gular size mitigates the impact of residual systematic uncertainties in modeling the effects of 
Suzaku's complex point-spread function (PSF), making the Perseus Cluster an ideal target to 
study cluster outskirts. 

A large mosaic of Suzaku observations of the Perseus Cluster, with a total exposure time 
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of 260 kiloseconds, was obtained in August/September 2009. The pointings extend along two 
arms from near the cluster center towards the east (E) and northwest (NW), out to a radius of 
2 degrees, which corresponds to 2.8 Mpc for a Hubble constant H = 70 km/s/Mpc. Here, we 
focus on the data obtained with the three available X-ray imaging spectrometer (XIS) cameras 
(see supporting online text and Fig. 1). We extracted spectra from annuli centered on the 
cluster center. After accounting for background emission, we modeled each spectral region as a 
single-temperature thermal plasma in collisional ionization equilibrium, with the temperature, 
metallicity, and spectrum normalization as free parameters. 

The best-fit radial profiles of temperature and metallicity are presented in Fig. 2. Individual 
elements are assumed to be present with Solar relative abundances (12). For comparison, we 
also show results previously obtained from an ultra-deep Chandra observation of the cluster 
center (13). The Suzaku and Chandra data sets show excellent agreement where they intersect, 
and together measure the temperature and metallicity structure of the intra-cluster gas with high 
precision and spatial resolution out to the virial radius (defined here as r 20 o, the radius within 
which the mean enclosed mass density of the cluster is 200 times the critical density of the 
Universe at the cluster redshift). In the narrow interval spanning 0.95-1. 05r 2 oo> the temperature 
is approximately a third of the peak temperature. Along the eastern arm, between 0.1 and 
0.7 Mpc, the temperature is systematically lower than towards the northwest, and the X-ray 
emission is brighter. This thermodynamic feature is known as a "cold front" and typically 
arises following a merger between the main cluster and a smaller sub-cluster (14). 

Our results show that the cluster outskirts are substantially metal-enriched, to a level amount- 
ing to approximately one third of the Solar metallicity. Previously, the only measurement of the 
metallicity close to the virial radius was obtained from a large region spanning 0.5-1 r 2 oo near 
the compressed outskirts of two interacting clusters (15). This previous result is in agreement 
with our measurements when converted to the Solar abundance units (12) adopted here. 
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From the Suzaku data, we have also determined the electron density, entropy, and pressure 
profiles, corrected for projection effects under the assumption of spherical symmetry (Fig. 3). 
Outside the cold front at 0.7 Mpc, there is a good match between the E and NW profiles, with the 
electron density decreasing steadily with radius, approximately following a power-law model 
n e oc r~ a with slope a = 1.68 ± 0.04. This is consistent with previous results from ROSAT 
data extending out to ~ 1.4 Mpc (16). 

Standard large scale structure formation models show that matter is shock heated as it falls 
into clusters under the pull of gravity. Simple theoretical models of this process predict that the 
entropy K should behave as a power-law with radius, K oc r /3 ', with /3 ~ 1.1 — 1.2 (17, 18). 
Except for the eastern cold front region, the entropy profile in Perseus roughly follows this 
expected trend until 2/3 r 200 . Beyond this radius and until 0.95r 20 o, both the E and NW arms 
show a flattening from the power-law shape, confirming hints from previous Suzaku results (7). 

The pressure profile is the most regular of the thermodynamic quantities plotted, and at most 
radii shows good agreement between the E and NW. At large radii, the pressure profile appears 
significantly shallower than expected by extrapolating the average profile of a sample of clusters 
studied previously with the XMM-Newton satellite within ~0.5r 200 (19). 

Invoking hydrostatic equilibrium of the ICM, the gas pressure can be used to estimate the 
underlying gravitational potential and total (dark plus luminous matter) mass profile of the 
cluster. Numerical simulations show that the latter typically follows a functional form described 
by Navarro, Frenk & White (20), also known as the NFW profile. We used the data from the 
northwestern arm of the Perseus Cluster, which appears dynamically relaxed, to determine the 
best-fit total mass profile, assuming an NFW form (supporting online text). The best-fit mass 
model parameters are typical of those predicted from numerical simulations; the NFW model 
provides a good description of the Suzaku data. 

Measuring the total mass profile allowed us to calculate the virial radius of the cluster, r 2 oo = 
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1.79 ± 0.04 Mpc, the corresponding enclosed total mass M 2 oo = 6.65_o 4 6 x 10 solar masses, 
and the cumulative gas mass and gas mass-to-total mass fraction, f gas , as a function of radius 
(Fig. 4). At relatively small radii of 0.2-0.3r 20 o, the measured f gas value is in good agreement 
with direct measurements from the Chandra X-ray Observatory (5) and measurements of the 
Sunyaev-Zel'dovich (SZ) effect (21) for two large samples of galaxy clusters. At about half 
of r 2 oo, the integrated gas mass fraction reaches the cosmic mean value computed from the 
CMB (22), considering that on average 12% of the baryons are in stars (23, 24) and the rest 
are in the hot X-ray emitting gas phase. Outside 2/3 of the virial radius, where the entropy 
also deviates from the expected power-law behavior, we find that the apparent f gas exceeds the 
cosmic mean baryon fraction measured from the CMB (22). 

The most plausible explanation for this apparent excess of baryons at large radius is gas 
clumping. In the X-rays, the directly measurable quantity from the intensity of the bremsstrahlung 
emission is the average of the square of the electron density, (n 2 e ), rather than (n e ). If the density 
is not uniform, i.e. the gas is clumpy, which is expected to occur as matter falls into the cluster, 
the average electron density estimated from the bremsstrahlung intensity will overestimate the 
true average, affecting the gas density, gas mass fraction, entropy, and pressure profiles. 

Outside the central region, and inside the radius where clumping becomes important, the 
measured f gas profile shows good agreement with recent numerical simulations (25), where a 
semi-analytic model was used to calculate the energy transferred to the intra-cluster gas by su- 
pernovae and AGN during the galaxy formation process. This model did not include a realistic 
implementation of gas cooling, and does not capture the complex processes in the central cool 
core of the cluster; the model is therefore not plotted in this region. Extrapolating this model 
into the outskirts where clumping is important, we used its predictions together with the mea- 
sured f gas to determine by how much the electron density must be overestimated to produce the 
difference between the data and the model. This factor (plotted in green in the bottom panel 
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of Fig. 4) reaches a value of up to 4 in the last annulus centered around the virial radius. The 
dense clumps are likely to be infalling and may be confined by ram pressure. 

Correcting the electron density using this factor, and accordingly the entropy and pressure 
profiles, we obtained the red lines shown in Fig. 3. The dumping-corrected entropy profile 
along the northwestern arm is consistent with the expected power-law profile. Moreover, the 
dumping-corrected pressure is also consistent with that expected by extrapolating the average 
profile of a sample of clusters previously studied with XMM-Newton (19). The corrected elec- 
tron density decreases more steeply with increasing radius, with the best-fit power-law index 
becoming a = 2.5. Correcting for clumping therefore seems to offer a simple solution to all the 
potential puzzles posed, at first glance, by the observed thermodynamic profiles of the Perseus 
Cluster at large radii. No additional physics is required by the data. 

Our study shows no evidence for the puzzling deficit of baryons at r > 0.5r 2 oo inferred from 
some previous studies of other systems, using lower quality data and/or extrapolated models 
(26-30). This suggests that within r < 0.5r 20 o the physics of the X-ray emitting gas is relatively 
simple and X-ray measurements can be used robustly for cosmological work. At larger radii, 
however, the cluster gas is significantly clumped. 

Numerical simulations predict gas clumping in the cluster outskirts (31). However, the 
amount of clumping depends on a large number of physical processes in the ICM which are 
currently uncertain, for example viscosity, conduction, star formation feedback, and magnetic 
fields. Although our results were obtained for just one galaxy cluster, it is expected that the 
observed physical processes are common. Our results therefore provide an anchor for numerical 
models of ICM physics and for simulations of the formation and ongoing growth of galaxy 
clusters. An independent measurement of gas clumping can be obtained from the combination 
of X-ray and SZ observations, which have different dependences on the electron density. 
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Figure 1 : X-ray surface brightness image of the NW (top) and E (bottom) arm mosaics observed with Suzaku, 
corrected for vignetting and instrumental background. The dashed white line marks the virial radius; the red circles 
mark excluded point sources and instrumental artefacts. The images have been rotated so that the cluster center is 
towards the left. 
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Figure 2: Projected temperature (kT) and metallicity (Z) profiles of the Perseus Cluster. Results from Suzaku ob- 
servations of the northwestern arm are shown in red, and of the eastern arm in blue. Earlier Chandra measurements 
of the cluster center {13) are shown in black. 
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Figure 3: Deprojected electron density (n e ), entropy (K), and pressure (P) profiles towards the northwest (red 
data points) and east (blue data points). The red line shows the northwestern profiles corrected for clumping. The 
expected entropy profile from simulations of gravitational collapse (17, 18) is a power-law with index j3 ~ 1.1, 
over-plotted as a black dotted line in the entropy panel. The average profile of a sample of clusters previously 
studied with the XMM-Newton satellite within ^0.5r 2 oo (19) is shown with a solid black curve in the pressure 
panel; its extrapolation to r2oo is shown with a dotted black line. 
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Figure 4: The integrated, enclosed gas mass fraction profile for the NW arm. The cosmic baryon fraction 
from WMAP7 (22) is indicated by the horizontal dashed black line; accounting for 12% of the baryons being 
in stars (23, 24) gives the expected fraction of baryons in the hot gas phase, shown as a solid black line. The 
values previously measured for a sample of relaxed clusters at smaller radii with Chandra (5) are shown with blue 
dots. Predictions from numerical simulations (25) are shown in green. The bottom panel shows by how much the 
electron density should be overestimated in each annulus due to clumping in order for the cumulative f gas not to 
exceed the correspondingly colored curves in the plot above. 
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Supporting text 

Observations and data analysis 

Below, we describe the observations and data analysis methods in more detail. The Suzaku 
mosaic of the Perseus Cluster consists of seven pointings extending towards the east and seven 
towards the northwest of the cluster center. The nominal exposure times moving outward along 
each arm are 10, 10, 10, 20, 20, 30, and 30 kiloseconds (ks). 

The standard screening criteria recommended by the XIS instrument tearrQ were applied 
to create a set of cleaned event lists. We removed times of low geomagnetic cut-off rigid- 
ity (COR< 6 GV) and verified that our observations were not affected by solar wind charge 
exchange (SWCX) emission using the method described in (SI). Images and spectra of the 
instrumental background were determined from dark Earth observations and subtracted from 
the data. The vignetting correction was performed using ray-tracing simulations of extended, 
spatially uniform emission (52). 

We extracted and co-added images in the 0.7-7 keV energy band for each pointing using 
the three XIS detectors. Spectra were extracted from annuli centered on the cluster center, 
excluding the visually identified bright point sources and instrumental artefacts marked by red 
circles in Fig. 1 of the manuscript. We also excluded a 30" region around the detector edges 
and around the calibration sources at the detector corners. 

The spectra were binned to a minimum of one count per bin and fit in XSPEC (S3) version 

12.6 using the extended C-statistic. We included in the spectral fitting data in the range 0.7-7 

keV and where the predicted contamination by stray light from the bright cluster center was 

below 20% of the best-fit cluster signal. For the Galactic absorption, we used the average 

column densities from the Leiden- Argentine-Bonn radio HI survey (S4). The cosmic X-ray 
1 Arida, M., XIS Data Analysis, http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/node9.html (2010) 
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background (CXB) was modeled with four components: three thermal models accounting for 
the local hot bubble, the Galactic halo, and a hot foreground present at low Galactic latitudes, 
and a power-law to account for the integrated emission of unresolved point sources. Data from 
the Suzaku mosaic beyond the estimated virial radius of Perseus, as well as from the ROSAT All 
Sky Survey (RASS), were used to constrain the parameters of the CXB model and estimate their 
spatial variations. Systematic uncertainties due to expected variations of the CXB are typically 
below the la statistical error in all the profiles shown in the manuscript. After accounting 
for background emission, each individual spectral region was modeled as single-temperature 
thermal plasma in collisional ionization equilibrium, using the apec code (55). 

NFW mass model 

The Perseus Cluster is a relatively relaxed system, as indicated by its regular X-ray morphol- 
ogy and moderately strong cool core. There is no evidence for an ongoing major merger either 
from X-ray data or from the internal galaxy distribution. A minor merger is readily identified in 
the thermodynamic profiles along the E arm, but no similar signatures are seen toward the NW. 
The NW arm is thus a particularly relaxed direction in a cluster that is relatively relaxed overall. 

To constrain the mass profile for the cluster, we used the Suzaku data for the NW arm and 
Chandra measurements for 4 inner annuli (0.02-0.08r 20 oX which capture the decrease in ICM 
temperature at small radii. Due to residual cross-calibration uncertainties (S6), the statistical 
precision of the Chandra measurements was limited not to exceed that of the Suzaku data, 
effectively using only a fraction of the available Chandra data. 

The mass analysis employs the standard assumptions of hydrostatic equilibrium and spher- 
ical symmetry (the latter is a good approximation when modeling the X-ray emission from 
relaxed galaxy clusters; 57). We parametrized the mass distribution using an NFW model, 
which is well motivated by simulations of cosmological structure formation (S8). The cluster 
atmosphere was modeled as piecewise isothermal (constant temperature in each shell). Apply - 

15 



ing the hydrostatic equation to the set of temperatures in each annulus, along with the NFW 
mass model, specifies the density profile of the ICM up to an overall normalization factor. This 
non-parametric description of the ICM allows considerably more freedom than methods which 
require the temperature and density profiles to follow simple analytic functions. This is particu- 
larly important when modeling the outskirts of clusters, for which no robust prior expectations 
for the forms of the gas density and temperature profiles exist. 

This model was implemented in XSPEC as an adaptation of the CLMASS mixing model 
introduced by (59). The predicted model emission from each spherical shell is projected and 
fitted to the spectra for all annuli simultaneously. The fitting procedure used a Markov Chain 
Monte-Carlo algorithm to determine posterior distributions for the gas mass, total mass and gas 
mass fraction profiles. The error bars shown in the upper panel of Fig. 4 are the marginalized 
(one-dimensional) 68.3 percent confidence intervals for the integrated f gas values at the outer 
edge of each isothermal shell. 

The minimum C-statistic (corresponding to the best-fitting parameters) is within the range 
expected from Monte Carlo simulations of the data, indicating that the model provides an ac- 
ceptable fit. The best-fitting NFW function has a concentration parameter c = 5.0 ± 0.5, which 
falls comfortably within the expected range for clusters of the measured mass, as determined by 
numerical simulations. The measured mass also agrees well with expectations based on simple 
mass -temperature scaling relations (e.g. S10). The gas temperatures and densities from the fit 
are in excellent agreement with the values determined from simple spectral deprojection. 

Our analysis shows that gas clumping becomes significant in the outer 5 annuli. Relaxing 
the assumption of hydrostatic equilibrium in these regions, while still accounting for the impact 
of projected emission appropriately, leads to consistent results on the gas mass fraction profile. 
Similarly, excluding the central 4 annuli (those using Chandra data) does not significantly affect 
the results; nor does relaxing the mass model to a generalized NFW form, which includes the 
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central density gradient as an additional free parameter. 
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